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ABSTRACT 

The hyperspherical harmonics (HH) provide a complete basis for the expansion of 
atomic wave functions, but even for two particles the number of harmonics for a given 
order is not trivial and, as the number of electrons increases, this degeneracy becomes 
absolutely prohibitive. We modify the method by selecting a subset of the basis that, 
we feel, will yield the physically most important part of the wave function, and test 
the idea on simple systems. 

In a previous work (with M. Haftel) of the ground and first excited 1 Sq states of 
the helium atom we found that the error in the binding energy of the 1 Sq ground state 
was of the order of 1 part in 10 4 and that it decreased for the first excited state. 

We now have applied our method to the equivalent triplet 3 Si states. We report 
on this work, and our results, and hope to draw attention to the interesting accuracy 
that we obtain with the relatively simple wave functions of our formulation. 



Introduction 



The bound states of the helium atom have been studied by using several methods. 

The most accurate results have been achieved by variational calculations where the 
states are described by wave functions which include more than 200 parameters 
Other, more recent, calculations have been performed using expansions in hyperspher- 
ical harmonics (HH), where the convergence was speeded up by the introduction of a 
correlation factor describing the average behaviour of the wave function at short and 
long ranges0. The residual part of the wave function was expanded in HH. This leads 
to an accuracy of about nine figures in the binding energy, but with a nonuniform 
convergence as the number of HH is increased. 

In our approach to the many body problem, we select the HH occurring in the wave 
function by following a scheme in which the basis elements are introduced according 
to the degree of complexity of the correlations that they describe. 

In the first approximation, a HH of minimal degree (because of the Pauli princi- 
ple) defines the state. In the next approximation, terms will appear in our equations 
consisting of the product of the HH of minimal degree and the potential. The basis el- 
ements that can correctly reproduce these terms form our Potential Harmonic basis 0, 
which we use in this work. A further approximation would be obtained by applying 
the potential to the wave function, including the correlations already obtained from 
the previous order and by taking into account the new basis elements generated by this 
product. 

In the analysis of the ground and first excited 1 5'o state of the helium atom (H, what 
was studied was the accuracy obtained by using the PH basis. The resulting expansion 
of the wave function provides a simple structure where one and two body contributions 
are exactly separated. It leads to a simple description of the wave function. 

We found that the error in our calculation of the binding energy of the 1 Sq ground 
state was of the order of 10~ 4 and decreased for the first excited state. 

Our purpose, then, in this paper was to study the worth of the method for triplet 
3 5i states. Since the potential basis is complete in so far as describing the independent 
motion of the electrons (i.e. when the interaction between the electrons is neglected), 
the missing energy must come from the Coulomb interaction between the two electrons. 
In the 1 S'o state this interaction is maximum at short distances and is expected to be 
larger in the space symmetric 1 S'o states than in the space anti-symmetric 3 Si states, 
where the effect of this interaction damps out at small distances where the wave function 
goes to zero. We were therefore very hopeful that our method would be even more 
satisfactory for the triplet than for the singlet. 
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Formalism 



In the HH expansion method the coordinates are expressed in polar form in a 
D-dimensional space where D = 3N is the number of independent space degrees of 
freedom. N = 2 for two electrons with position vectors x± and x*2, relative to the nucleus 
with nuclear charge Z and polar coordinates (rj,a;j), % = 1, 2, where the cJj denote the 
two standard angle measures for the ith electron, the hyperspherical coordinates (r, Q) 



are chosen to be the hyperradius r = yrf + r^ and the set Q of five angular coordinates, 
uji, uj 2 and z, where r\ = r cos0, r 2 = r sin0 and we replace <ft by 



COS. 



r\ — r\ 



The Schrodinger equation for the helium atom, in atomic units (h = m = e = 1), 

|-~(V? + V 2 2 ) - Z(- + -) + J- - e\* = , (2) 
12 ri r 2 r i2 J 



n r 2 r 12 
is then transformed, in hyperspherical coordinates, into 



d 2 L 2 {9) 




Vl^y r Jl-xy/l 



z 2 



E}$(r,n)=0 (3) 



where a; is the cosine of the angle between the two electrons, (x = xi.x 2 /(rir 2 ) = cos 9), 
L 2 (Q) is the grand orbital operator, and \l/ and $ are related by 

^ = r - 5 / 2 $(r,fi) . (4) 

We note that if Y[n(f2) is a HH, with r L Y[£](f2) the associated harmonic polynomial of 
degree L, then L 2 (Q) in equation (|3|) satisfies 

[L 2 (Q) + C(C + 1)]Y [L] (Q) =0 (5) 

with £ = L + 3/2in the six dimensional space spanned by the coordinates x\ and x 2 . 
Also, the Hamiltonian in (Q) is an even function of z and therefore the wave function 
can have a definite z-parity. The ^-parity is associated, according to ([I]), with the 
parity in the exchange of the two electrons, i.e. it is even for the singlet 1 S$ states and 
odd for the triplet 3 S\ states. 

When the interaction between the electrons is set to zero, equation (|3|) is separable 
and the wave function for S'-states can be written either as a product of individual 
wave functions 

* = <Pmh (n) <Pn 2 l 2 (r 2 ) ± <J>n 2 l2 ( r i) <Pmh ( r 2) (6) 
where +(— ) is for space (anti) symmetric wave functions, or as a function 

with a z-parity even(odd) for space (anti) symmetric wave functions. Both descriptions 
are equivalent. 
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Let the HH expansion of the wave function in ([|) be 

$(r,n) = J2Y [L] (n)u [L] (r) 



where [L] is the set of the five quantum numbers [L, li, l 2 , mi, 1712] associated with the 
five degrees of freedom (uji,uj 2 ,z). The (h,mi) are the orbital quantum numbers and 
L = 2k + li + l 2 where k is the quantum number associated with z. 
The coupled equations for the partial waves U[L](r) are 



d 2 C(C + 1) 



+ 



E ) «[L]( r ) + E < Y lL]\V(r,Q)\Y [L/] > u [v] (r) = (9) 

J [L'\ 



where in our problem V(r, Q) is the potential occurring in equation (|3p. 

In equation @ the sum starts from a minimal degree L m . In the hypercentral 
approximation (HCA) only the first term in equation (|8]) is taken into account and 
equation @ is reduced to a hyperradial equation 



i - 

[2 dr 

where 

For a normalized HH 



+ U [Lm] (r)-E\ U[Lm] (r) = 0. 



U [Lm] (r) =<Y [Lm] \V(r,n)\Y [Lm] > . 



(10) 



dtt 



Here, the surface element of the unit hypersphere is 

dVt = du\ dtOi (sin cos < 



duj\ duj 2 - Vl — z 2 dz . 



For S'-states 
and 



[2fc,0,0,0,0] 



(Q)=n-^ 2 Cl(z) 



(12) 



dtt = tt 2 sin 9 d6 Vl-z 2 dz . 



With the Gegenbauer polynomials Cq = 1 and C\(z) = 2z, the hypercentral potential 
for the 1 5'o and 3 Si ground states are, for k — 0, 1 respectively, 



and 



- I V(r, Q) sin9d6Vl-z 2 dz 

7T J 

- / V(r,Q) sm9d9z 2 Vl-z 2 dz 
r J 



(13) 
(14) 



where the integrals are taken in the domain < 9 < n and —1 < z < 1. 
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For a Coulomb potential, the integrals fll^D and are simply proportional to 
1/r, and therefore the HCA equation ( |ToD can be solved analytically ||. In order to go 
further, another variable must be introduced. An expansion of in terms of Legendre 
polynomials is obtained from the generating function 



Y J Pi{u)/ = {l-2zu + z 2 )- 1 '\ (15) 



leading to 



1 



r 12 r y/l — sin 20 cos# 

We expand the wave function for electrons in ^-states in terms of normalized Legendre 
polynomials: 

oo 

*(r, z, x) = J2?e( x ) u t( z > r ) , ( 17 ) 



where P(x) = J £ + l/2?<(x). Introducing this form for the wave function into Eq. 
(0), we obtain an infinite system of coupled second order partial differential equations 
for the partial waves ue(z,r): 



d 2 L%z) 

q^2 + -^2— + 2 ^ 2 - 



1 

+ 



+ 2E\ u £ (z,r) 



OO 

2^ <P e \ — \P e >u e {z,r) (18) 



£'=0 



Tl2 



where 



and the matrix element 



where is given by (fl6|) in terms of r, 2 and x = cos 9. Rather than try to solve this 
system of coupled two-variable second order partial differential equations to atomic 
accuracy, we expand $(r, z, x) in hyperspherical harmonics. To then solve exactly 
equations fll8|) requires the complete Zernike-Brinkman basis || HE for two particles 
in S-states, i.e. a basis complete in the variable z for each value of i in the expansion 
(|T7l). This would lead us to a large number of coupled equations (|9]) to be solved for 
a given accuracy and we try to reduce the size of the basis and choose the potential 
basis. 
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This basis consists of: (i) the basis for t = which is complete for an expansion 
of any function of z only (in particular of the potential rj" 1 + r^ 1 ), and (ii) the basis 
needed for a complete expansion of r^ 2 . This implies that equation ([|) can then be 
solved with a harmonic expansion which includes only two orthogonal elements for each 
grand orbital L = 2k, except when L = 0. 

The PH approximation selects the partial waves in expansion (H) which are coupled 
to the largest components, L — (L — 2), of the 1 S ( 3 «Si) ground states of the wave 
function respectively We call this the optimal subset. 

In order to extract this subset, we introduce the kinematic rotation vector R(tp) = 
fi cosy? + r 2 siny?, with parameter ip, where -R(O) = ri, R(tt/2) = r 2 , .R(37r/4) = 
(f 2 — ri)/V2, and we define z(cp) = 2R 2 ((p)/r 2 — 1, with z(0) = cos 20 = z, leading to 
z((p) = zcos2(p + xy/l — z 2 sin2<£> and 

1 y/2 1 V2 1 1 



VI + ^ ' r 2 VI -z' r i2 r ^/l + z(3tt/4) 



(20) 



for = 0, 7r/2, and 37r/4, respectively 

The Potential Harmonics for the S'-state are given in terms of the kinematic angle 

K(^V) = ^ /2 Cl{z^)) (21) 

where C\ is a Gegenbauer polynomial of even grand orbital L = 2k. We use the notation 
"P 2fc (f2i), i — 1) 2, 12, for P 2fc (f2, </?j), where <£> = 0, 7r/2, 37r/4 refer to the respective 
choices of the kinetic angle parameter for each coordinate. The PH expansion of rj , 
for i — 1, 2, 12, is then given in |3j for n = — 1 by 

n r V (2fc + l)(2fc + 3) 

where V^i^) = ( — l) fc ^2A:(^i)> according to the parity of the Gegenbauer polynomial, 
and 

1 _ 8V2T v fe fc + 1 Q 
r^ - — V ] (2k + l)(2k + 3) V ^ Ql2) ■ 

In a PH expansion of the wave function for the singlet S-state, the optimal subset 
consists of: (i) V$ k (£li) for k even and (ii) that part of V^i^u) orthogonal to V^iP-i)- 
One uses the orthonormalized base B k {VL) = V% k (fli) for k even, B k (VL) = V^i^n) for 
k odd, and B£(Sl) = ^7== [^(^12) - (^1)^(^1)] for k even > where the matrix 
elements 

<n k {n 2 )\v" 2k {Q 1 ) >= (-1 

and 

<^2 fe (^)l^(^12)> = 



are used in agreement with equation (25) of reference 




k odd 
k even 
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For triplet 5-states the B k (Q) for k odd are suitable for the motion of the two 
electrons around the nucleus, but V 2k (Qi 2 ) is not appropriate since it is even in the 
exchange of the two electrons. In order to find the second set of HH we note that, since 
the normalized PH are proportional to Gegenbauer polynomials (0), then 

<V^ k {Q,f)\V^,{Q,f) >=5 k)k , . 

By taking the derivative with respect to (p, noting 

one generates the additional basis elements 

dv° 2k (n,<p) 2 



Of 7T 3 / 2 



Vl — z 2 cos 2p — z sin 2p C\_ x (z cos 2<p + x\J\ — z 2 sin 2pj 



For if = one obtains the set of HH where electrons in a p-state are coupled to 
form an Estate obviously orthogonal to B k (Q), where both electrons are in an s-state. 

For if = 37r/4, one obtains the required basis zC\[x\/\ — z 2 ), which however must 
be normalized and has to be made orthogonal to B k (Q) for k odd. The desired basis 
elements are 



^—zCl^x^fY^z 2 ) fc>0even 



7T 3 / 2 y/k{k+2) 
2^3 



TT 3 / 2 y/k(k + 2)-3 



k > 1 odd 



It is worthwhile to note that, since z(37r/4) = r\ 2 /r 2 — 1, the basis element B k 
describes the electron pair in an S'-state. Then, in the expansion 



k,k i 



the electrons are individually in s-states in B k . (Here and in the appendix, we some- 
times emphasize that indices take on the values that are associated with the perpen- 
dicular set by using the notation k± or k' ± . Where there is no room for confusion, we 
drop the _L.) This means that in our approximation the space wave function becomes 
the sum of two functions. In the singlet state 

= ips(rl - rf, r) + 4>j(r 2 2 , r) 

where ip$ and ipg are even in an exchange of the two electrons for the singlet state, 
while in the triplet state 

#t = ^ T (r\ - rl, r) + [r\ - r\)^{r\ 2l r) 

where tpx is odd in the exchange of the two electrons. 
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The partial waves are solutions of the system of coupled second order differential 
equations 

(T k - E)u k + ]T {Utuv + U^u^) = (22) 



k' 

i {U k u k , + C 

A-'./,- 

and 



(T fc± - + £ + Uiu w ) = (23) 

where 



n = _l_ + £ fc (£ fc + f) 
o?r 2 r 2 
and C k — 2/c + 3/2, (both valid for and k±). 

The matrix elements of the potential matrix are given in Appendix A. 
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Results and Analysis 



The r dependence of the coupling matrix is simply 1/r (i.e. multiplying a matrix 
of constant coefficients). To solve the resulting coupled differential equations, we take 
advantage of a method developed by V. Mandelzweig and M. Haftel[|| - and embodied 
in one of their numerical programs. It was used by Haftel and us (here, Fabre and 
Larsen) to obtain results for the singlet [[J]. Michael kindly made available to us his 
program, that we used for the triplet. 

Its essence is to calculate the wavefunction and its logarithmic derivative by using 
expansions - and, then, obtain the matrices of coefficients (of these expansions) by 
numerical recursions. 

We solve the set of equations for a number of values of k max , such as 33,37,41 ... 
and extrapolate the result to infinite order. We can compare our results to a wealth of 
reliable data|], |J. 

Singlet 

As noted above, the results have been published |3[]. The result for the ground state is 
that we obtained an extrapolated value of 2.903471 atomic units, compared to the more 
accurate variational result of 2.903724. (We truncate all our numbers to 6 decimals.) 
The error was thus ~ .01%. 

For the excited state, the extrapolated value was 2.145825 compared with a varia- 
tional value of 2.145974. The resulting error was thus ~ .007%. 

These were based on calculations for k max = 24, 28, 32, 36 and noting0 that, asymp- 
totically, the differences between the energies for successively larger values of k max de- 
crease at least as fast as k^ ax . The differences were then fitted to an inverse polynomial: 

AE(k) = E(k) — E{k — 4) 



k 2 (A + Bk + Ck 2 ) 



and the extrapolated values were obtained by adding the remaining contributions that 
would result from summing to infinity. 

Another point has, in retrospect, proved important. For the ground state, the 
value for k max = 36 is 2.903067, i.e. this 'raw' results is already correct to 4 digits. The 
extrapolated 'tail' only adds a modest amount. For the excited state, the situation is 
quite different. The 'raw' value is 2.129457 and we see that therefore the difference 
between the resulting energy and the correct value is, roughly, 1%. The contribution 
of the extrapolated 'tail' is very important yielding several more digits in agreement. 



Triplet 

In Table 1, we present the numerical results (in a.u.) obtained by solving the coupled 
equations. These need, again, to be extrapolated, so as to take into account the 
contribution of the higher order potential harmonics, and then compared with good 
variational estimates: 2.1752495 for the lowest state, 2.068685 for the excited state. 
Without the r 12 interaction, the exact result is 2.5 a.u. 
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Table 1: energies from solving the equations up to k. 



kmax 


ground state 


excited state 


no ri2 interaction 


21 


2. 


1401236235411021 


1.9201801253586101 




23 


2. 


1459023327465288 


1.9377802892305477 




25 


2. 


1504882412407837 


1.9525548534156408 




27 


2. 


1541746943943439 


1.9650792663836013 




29 


2. 


1571714798190449 


1.9757878236977571 




31 


2. 


1596321992027180 


1.9850138971543965 




33 


2. 


1616708453178672 


1.9930169776811079 




35 


2. 


1633735780590255 


2.0000017724075086 




37 


2. 


1648062043495217 


2.0061315935556941 




39 


2. 


1660197311520152 


2.0115381473492293 




41 


2. 


1670540234821436 


2.0163286338520296 




43 


2. 


1679406196279248 


2.0205910881663552 




45 


2. 


1687046390433931 


2.0243983641445255 




47 


2. 


1693662957017716 


2.0278112036995569 


2.4978184400838062 


49 


2. 


1699419469276808 


2.0308805802207102 


2.4980527934858051 


51 


2. 


1704449472143039 


2.0336495409553854 


2.4982549671085871 


53 


2. 


1708862505879244 


2.0361546418675421 


2.4984302656593832 


55 


9 


1 71 27491 227421 1 1 


2.0384270954208123 


2.4985829842420010 


^7 

O 1 


2. 


1716184476775823 


2 n4D4Q^fi7QS^7fiQ4Q 


9 4Q871 fifi9441 ^9fi1 fi 


59 


2. 


1719231375119347 


2.0423774774219713 


2.4988340592267531 


61 


2. 


1721942486212956 


2.0440984689458583 


2.4989376603658377 


63 






2.0456740179909737 


2.4990293968470676 


65 






2.0471192766548691 


2.4991109122218091 


67 






2.0484475085284343 


2.4991835853845747 


69 








2.4992485787876214 



The first thing that we notice, then, is that, even for the much higher values of k max 
that we have used for the triplet, the difference between our results and the variational 
ones is still more than 1 part in a thousand for the ground state and 1% for the higher 
state. If we had limited ourselves to the orders used in the singlet, our errors would 
have been 1/2% and 3%, respectively! 

We were thus driven to higher values of k max and thus to a greater number of states 
in the calculations. For example, a k max of 61 corresponds to a total number of states 
of 91. Even then, we also had to work harder at the extrapolation than we did for the 
singlet. The differences were fitted to an inverse polynomial with more coefficients: 



w Ak 4 + Bk s + Ck 2 + Dk + E 

with the As variously defined as AE(k) = E(k)-E(k-A) and AE(k) = E(k)-E(k-2) 
and the extrapolated values obtained by using this fit, for the differences, to obtain the 
higher order contributions. 
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Examination of the results of the extrapolations then shows that the best strategy 
consists in fitting the coefficients with the minimum of data points, taken at the highest 
values of k max . For the ground state, the use of A(k) = 4 yields the best result. For 
the excited state, A(k) = 2, is superior. 

Thus we find that for the ground state, the use of the set of values of k max \ 
(51,53,55,57,59,61) yields 2.175042, while (41,45,49,53,57,61) gives us our best result 
2.175144 (both results, here, rounded to 6 decimal places). We recall that this must 
be compared with 2.1752495. 

For the excited case, (55,57,59,61,63,67) gives 2.068704, while (47,51,55,59,63,67) 
gives 2.068757, also rounded at the same place. We compare this with 2.068685. 

Since in our method we only approximate by truncating the basis, both the elements 
which appear in it and the highest order that we take into account in our calculations, 
our results must obey the variational principle and must yield upper bounds in the 
(negative) energies of the states. We can verify that this is indeed the case for our 
"raw" results, displayed in Table 1. For the extrapolated results, this is only evident 
for the ground state. Our results for the excited state are slightly deeper than the 
correct energy. 

This is an artifact of the extrapolation. Clearly, our inverse polynomial only im- 
perfectly represents the correct asymptotic behaviour of the differences. In fact taking 
different sets of points, from which to extrapolate, and trying to 'see' trend lines, we 
think that we see oscillations and certainly instability and serious limitations in our 
ability to extrapolate. 

Nevertheless, we obtain energies to within 1 or 2 parts in 20, 000 for the lowest state 
and 1 to a few parts in 100, 000 for the excited state. 

Finally, we note an extrapolation with A(k) of 2, for the case where the electrons are 
uncoupled. Our result is then 2.499999995, rounded to nine decimal places, compared 
to the exact value of 2.5 In this case our basis is complete, aside from truncation in 
the order, and we take this result as a confirmation of the soundness of our numerical 
approach and procedures. 
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Conclusion 



We believe that we have demonstrated that the use of Potential Harmonics is a 
sound way of obtaining reliable results, at least for 2-electron atoms, to an accuracy 
of 1 part in 10 4 or 10 5 . Our aim was to confirm this, with our first two calculations, 
and then to proceed with higher excited states for these 2-electron atoms and, perhaps 
more excitingly, with 3-electron atoms! 

While this is probably feasible, and in fact the potential harmonics seem to yield 
better accuracy for the excited states as 3-body and higher correlations become less 
important, our experience leads us to the following caveat. 

We have found many aspects of our effort, such as the calculation of our matrix 
elements, the necessary extrapolations, to be quite delicate! We would prefer other 
methods which, using integral equations, would not require them. 

Our final and best thought, then, is to draw attention to the simplicity of the form 
of the wave functions that we have used, and to the fact that, in them, the pairs are, 
each, in an s-state. 
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Appendix A 



The matrix elements in equations (|22[) and (|23|) , after making the potential harmonics 
basis substitution, take the form of integrals of products of three Gegenbauer polyno- 
mials. Let us introduce a number of expressions for sums that occur from particular 
combinations of HHs and PHs: 



k+k' 

s k k = E 

m=\k— k'\ 
(step 2) 



(m + 1) 



(2m + l)(2m + 3) 



where Si = S?' 



Or 



k+k' 

E 

m=\k— k'\ 
(step 2) 



(2m + l)(2m + 3) 



and o\, = erf. 



E 



y 



^A^ 1 Jk + k' -m + l)(k-k' + m + l) 

> (m+1)- ; 77 r -x 

^ K ' (2m + 5)(2m + 3) 

m=[fc-Jfc'|+l 
(step 2) 

(k' - k + m + l){k + k' + m + 3) 
(2m + l)(2m- 1) ' 



with £*, = E£'. 

(E x )r 



x 



16 ^ (_1)n2 (2n + l)(2n + 3) ^ 

n=\k-k'\ m=\k-k'\+l 
(step 2) (step 2) 

(m - fc + k' + l)(m + fe - fc' + l)(/c + fe' + 1 - m)(/c + k' + 3 + m) 

m(m + 2) 

where, by symmetry, (S- 1 )^, = (E- 1 )^'. 
The matrix elements 

First, consider the contribution from the two electrons separately around the nucleus, 
where k and k' are both odd because of the anti-symmetric requirement of the triplet 
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wave function. 



Tjk> _ 

u k — 



r 

64 

7rr 



Bk 







Z 


z , 


— Bk 1 ') 


- 2 -( Bk 




+ — 


ri2 1 


r \ 


n 


r-2 



B^ 



4 



— a\ - ZS k 



k' 



(24) 



Then, consider the interactions between the electron-electron wave function and the 
electron-nucleus wavefunctions. Here, we can have k odd and k' ^ even, 



T k' 





2 /„ , Z 


Z , 


- (Bk — B^} 


>--(B k \- 


+ — 


r \ r 12 1 


r \ T\ 


r 2 



'^k ' 



nr^/k'(k' + 2) k(k + 2)' 

and for the case k' odd and k ^ even, we use = U k± . 
For the case where k and k' are both odd, 



(25) 



y 



16^/6 



irrJk'(k' + 2) - 3 [ k(k + 2) 



( " 1)( " 1)/2 eI' - (-D (fc " i)/2 -r 



(26) 



Finally, we need to consider the perp-perp matrix elements. Thus, we have for k and 
k' both even: 



u k '- - -Ib^\ — 



B t)-~-{B k 



- + -\Bt> 

n r 2 



16 



k and k! both odd, 



nr^k(k + 2) ^Jk'(k' + 2) 



16 



(27) 



|V2 

- 12 Z 

odd and /c' even, 



■Kr\jk(k + 2) - 3 ^k'(k' + 2) - 3 
3 3 



x 



\ k{k + 2) k'{k' + 2)) 

E ± )r-(-i) (w/2 ^']}, 



e:'+3(-i)i^'i/ 2 4 



(28) 



16\/2 



Ttr^k(k + 2)^Jk'(k' + 2) -3 
and finally, with k even and A;' odd, 



k'(k' + 2) 



16V2 



nryjk{k + 2)^/k'(k' + 2) -3 V + 2 ) , 



Eh! 
k ' 



Ek 
k! ' 



(29) 



(30) 
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